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PROPAGATION OF GAMMA-RAYS IN MASSIVE BINARY CEN X-3: INTERACTION OF 
CASCADE GAMMA-RAYS WITH THE MASSIVE STAR 
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ABSTRACT 



We consider the propagation of very high energy 
(VHE) gamma-rays in the radiation field of the mas- 
sive star of binary system Cen X-3. VHE gamma- 
rays or electrons, injected by the compact object 
close to the surface of a massive companion, develop 
inverse Compton e ± pair (ICS) cascades. Based on 
the Monte Carlo simulations, we obtain the fraction 
of secondary cascade photons which collide with the 
surface of the massive companion in Cen X-3 system. 
The distribution of photons falling on the surface of 
the massive star is investigated. These photons in- 
teract with the atmosphere of the star and should 
excite different 7-ray lines (nuclear, annihilation 
line). We estimate that these 7-ray lines can be likely 
detected by the INTEGRAL telescopes. 
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ries: Cen X-3. 



1. INTRODUCTION 



Observations of some massive binaries suggest that 
they may be sources of 7-rays in the GeV and TeV 
energy range. For example the massive X-ray binary 
system, Cen X-3, containing a neutron star with a 
4.8 s period in a 2.09 day orbit around an O-type 
supergiant, has been detected above 100 MeV 7-rays 
by the EGRET detector on the Compton Observa- 
tory (Vestrand et al. 1997). There are evidences that 
observed 7-ray emission at these energies has a form 
of outbursts and is modulated with a 4.8 s period of 
the pulsar. Based on the observations in late 80's, 
two groups (Brazier et al. 1990, North et al. 1990) 
reported Cherenkov detection of positive signal at 
TeV energies from Cen X-3 at an orbital phase of 
~ 0.75, which is modulated with a period of the pul- 
sar. This emission has been localized by Rauben- 
hcimer & Smit (1997) to a relatively small region 
between the pulsar orbit and the surface of a massive 
companion which may be the accretion wake or the 
limb of the star. More recently the Durham group 



has detected a persistent flux of 7-rays above 400 
GeV on a lower level than previous reports (Chad- 
wick et al. 1998, 1999a, b). No evidence of correlation 
with the pulsar or orbital periods has been found 
and no evidence of correlation with the X-ray flux 
has been detected (Chadwick et al. 1999a, b). 

The problem arises if such TeV photons injected 
close to the surface of the massive star (e.g. by 
the compact object) can escape from the soft ra- 
diation field of the massive star in the binary sys- 
tem Cen X-3. The computations of the optical 
depth for TeV photons in the radiation fields of 
such type massive stars show that the absorption of 
TeV photons should be significant (e.g. Protheroe 
& Stanev 1987, Moskalenko et al. 1993, Bednarek 
1997). Recently we have performed Monte Carlo 
simulations of cascades initiated by 7-ray photons 
and electrons with different energy and angular dis- 
tributions, which are injected by the discrete source 
in the radiation field characteristic for the Cen X- 
3 system (Bednarek 2000). The spectra of escaping 
7-ray photons and their light curves have been ob- 
tained and discussed in the context of recent GeV 
and TeV observations. In this paper we compute the 
7-ray power and spectra of photons which fall onto 
the surface of the massive star in Cen X-3. These 
photons should excite different 7-ray nuclear lines 
and also annihilation line. Their detection by 
INTEGRAL telescopes should give independent in- 
formation on high energy processes which seems to 
occur in Cen X-3 system. The following parameters 
for the massive binary Cen X-3 are assumed: the ra- 
dius of the O star r s = 8.6 x 10 11 cm, its surface 
temperature T s = 3 x 10 K, the binary star separa- 
tion a = 1.2 x 10 12 cm (Krzemihski 1974). 



2. GAMMA-RAYS FROM ICS CASCADE 



Let us assume that 7-ray photons are injected by 
the compact object which is on an orbit around the 
star with the parameters characteristic for Cen X-3 
system. These photons may create e pair in colli- 
sion with the soft star photon (Bednarek 2000). The 
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Figure 1. Schematic picture of the cascade initiated 
by a primary gamma-ray fy) in the radiation field 
of a massive star. Gamma-ray photon, injected by 
the compact object or produced by primary electron, 
creates e ± pairs (e) in the interaction with a soft star 
photon e. The secondary e ± pairs, by scattering soft 
star photons, create secondary gamma-rays (71, 72, 
73/ Some of them (js) may fall onto the surface of 
the massive star. 

secondary pairs can next produce ICS 7-rays, initi- 
ating in this way the ICS cascade in the massive star 
radiation which is seen anisotropic in respect to the 
location of the injection place of the primary photons 
or electrons. 

The propagation of secondary pairs inside the binary 
system is determined by the strength and the geom- 
etry of the magnetic field which in general may be 
very complicated. We analize the case in which the 
random component of the magnetic field isotropize 
secondary cascade e^ pairs close to the place of their 
creation before they next scatter soft photon from 
the massive star. This condition is met if the mean 
free path for ICS process is longer than the Larmor 
radius of secondary pairs in the local magnetic field 
(see discussion and condition for that to occur in 
Bednarek 1997). We neglect also the synchrotron 
losses of secondary pairs in respect to their ICS losses 
on star photons. This can be done for the values of 
the magnetic field below certain limit (see Eq. 3 in 
Bednarek 1997). The bremsstrahlung losses of pairs 
in matter with the star wind densities characteristic 
for the OB stars are negligible. 

In the previous paper (Bednarek 2000), we show the 
details of the Monte Carlo simulations and concen- 
trate on the photons which escape from the binary 
system. However significant part of secondary cas- 
cade photons should collide with the massive star. 
The distribution of these photons on the star sur- 
face, their spectra and the fraction of photon power 
which fall onto the surface in respect to the power of 
all produced cascade photons is investigated in this 
paper. The considered picture is shown in Fig. [l]. 

We concentrate on two cases which were already dis- 
cussed in Bednarek (2000), i.e. the case of isotropic 



injection of electrons and isotropic injection of pho- 
tons with a power law spectrum and index -2. Elec- 
trons with such distribution can be accelerated at 
the shock front created in collision of the pulsar 
wind with the surrounding matter as proposed in the 
model by Kennel & Coroniti (1984). If the compact 
object (neutron star) is surrounded by dense cocoon 
created by matter accreting spherically-symmetric or 
by the accretion disk, then the 7-ray photons should 
be injected into the volume of the binary system. 
Based on our simulations we show in Figs. |^a and b 
the distribution of secondary photons on the surface 
of the massive star and the distribution of the 7-ray 
power of these secondary photons as a function of the 
cosine angle a, measured from the direction defined 
by the centers of the massive star and the compact 
object. Note that this distribution is axially sym- 
metric (see Fig. 1). 

The numbers of secondary photons and their power 
falling onto the massive star drop strongly with the 
angle a. Only a small number of photons is able 
to fall onto the opposite side of the massive star 
than the location of the compact object (the source 
of primary particles initiating the cascade), i.e. at 
cosine angles less than cos a = 0. This angular dis- 
tribution of falling photons is much sharper than 
the angular distribution of cascade photons escap- 
ing from the system. These interesting features are 
due to the type of the considered cascade process 
in which the secondary pairs are isotropized by 
the random magnetic field. These pairs produce sec- 
ondary ICS photons moving towards the massive star 
surface with higher probability. However for sec- 
ondary photons it is easier to escape from the system 
on the oposite side of the massive star than collide 
with the surface of the massive star on the oppo- 
site side. There are also differences in distributions 
of secondary photons on the massive star surface in 
the case of injection of primary photons and elec- 
trons (compare dotted and full histograms in Figs. ||a 
and b). The distribution of secondary photons from 
cascades initiated by primary photons shows higher 
number of secondary photons falling onto the op- 
posite side of the massive star and on the surface 
of massive star just below the compact object the 
than in the case of injection of primary electrons. 
This is due to the fact that primary photons pro- 
duce first generation of pairs in quite extended 
region around the compact object in comparison to 
the point like injection of primary electrons in the 
second case. We show also in Figs. ||a and b the 
spectra of secondary photons falling onto the mas- 
sive star in both discussed cases for different range 
of cosine angles a. As expected, the secondary pho- 
tons with the highest energies fall mainly onto the 
massive star for small angles a. On the opposite 
side of the massive star fall only secondary photons 
with relatively low energies. 

We estimate also the ratio of power of secondary pho- 
tons which fall onto the massive star to power of 
photons which escape from the system. In the case 
of injection of primary photons this ratio is equal to 
~ 13% and in the case of injection of primary elec- 
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Figure 2. The number of ~f-rays (a) and 7 -ray power (b) falling on the surface of the massive star in Cen 
X-3 system. It is assumed that the point like source injects isotropically primary j-rays (dotted histogram) or 
primary electrons (full) with a power law spectra and spectral index -2. 
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Figure 3. The spectra of j-ray photons falling onto the surface of the massive star in Cen X-3 system within 
the cosine angle a: cosa = 0.9 «-» 1 (full histogram), 0.5 <-*■ 0.6 (dotted), 0. <-> 0.1 (dashed), and —0.5 «-> —0.4 
(dot-dashed). The primary electrons (a) and photons (b), with a power law spectrum and spectral index -2, are 
injected isotropically by the compact object at the distance of 1.4 radii from the massive star. 



trons it is equal to ~ 40%. In the case of primary 
electrons, this is much more than expected from 
simple geometrical considerations (the case of non- 
cascading photons) which predicts that only ~ 13% 
of photons should fall into the star surface since this 
part of the sky is shaded by the massive star. In the 
case of primary photons, the non-cascading photons 
dominate and the ratios are comparable. 

This results show that significant part of the 7- 
ray power observed from the binary system should 
be transferred onto the surface of the massive star. 
These 7-rays can heat the star and partially can be 
re-emitted in the form of nuclear lines and e^ anni- 
hilation line. Therefore we suggest that the binary 
systems of the Cen X-3 type can also produce 7-ray 
lines in the MeV energy range. 

To find out what is the relation between this ex- 



pected 7-ray line emission and the continuous emis- 
sion, we compare in Figs. f|a,b and c the spectra of 
photons which fall onto the surface of the massive 
star at fixed angles a with the spectra of secondary 
cascade photons which escape from the system at 
these same ranges of angles (see Bednarek 2000). 
Fig. ^a shows that the line emission should be the 
strongest when the compact object is in front of the 
massive star. Similar behaviour is expected for the 
escaping cascade 7-rays with energies above 100 GeV 
(ground Cherenkov telescopes) . On the contrary, the 
highest fluxes of escaping cascade photons with ener- 
gies above 100 MeV (space telescopes) are expected 
in the direction tangent to the massive star limb (see 
Fig. [|c) . Therefore it is expected that the maximum 
of the 7-ray line light curve should not correspond 
to the maximum of the 100 MeV 7-ray continuous 
emission. 
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Figure 4- Comparison of the secondary ~f-ray spectra which escape from the binary system Cen X-3 (full his- 
tograms) and which fall onto the surface of the massive star (dotted histograms) for these same cosine angles a 
equal to cosa = 0.9 <-> 1. (a), 0.5 <-> 0.6 (b), and —0.5 «-> —0.4 (c). The primary electrons with a power law 
spectrum and spectral index -2. are injected by the compact object. 



3. CONCLUSION 



Assuming that VHE 7-ray photons or electrons are 
injected inside the binary system Cen X-3, as postu- 
lated by observations of this source in the GeV and 
TeV energy ranges, we consider the cascades initi- 
ated by these particles in the radiation field of the 
massive star in Cen X-3. The spectra of secondary 
photons, which escape from the system and fall onto 
the surface of the massive star, are obtained by us- 
ing the Monte Carlo method. The photons which fall 
onto the massive star should excite the 7-ray nuclear 
lines and the e^ 1 annihilation line. The maximum in- 
tensities of 7-ray line emission is expected when the 
compact object is in front of the massive star. On 
the contrary, the maximum emission of the 100 MeV 
emission is expected for direction tangent to the limb 
of the massive star. Therefore the 7-ray light curves 
observed at low and high energies should appear dif- 
ferent. 

Our model predicts that about 13 <-> 40% of 7- 
ray power escaping from the system can be trans- 
ferred to the surface of the massive star as a re- 
sult of cascades in the binary system Cen X-3. The 
phase averaged 7-ray luminosity above 100 MeV, 
reported by the EGRET telescope, is ~ 5 x 10 36 
erg s _1 (Vestrand et al. 1997). Therefore about 
~ 1 <-> 2 x 10 36 erg s _1 should fall onto the star. 
The INTEGRAL Spectrometer instrument will have 
narrow-line sensitivity of a 5 x 10~ 6 ph cm~ 2 s _1 
at 1 MeV and continu um sensitivity 1.5 x 10~ 7 ph 
cm " 2 s -1 at 1 MeV (|http:/ /astro. estcc.esa.nl/SA~ 1 



gen |3ral/projects/iNTEGRAL| ). It will be able to dc 
tect the 1 MeV 7-ray source emitting the power 
~ 10 35 erg s _1 at the distance of Cen X-3 system 
equal to 8 kpc (Krzemihski 1974). Therefore the 
conversion efficiency of cascade photons into the 7- 
ray line photons on the surface of the massive star 
should be ~ 3 <-> 10%. This value does not seems to 



be extraordinary. The detailed predictions for pos- 
sible detection of such lines from Cen X-3 by the 
INTEGRAL telescopes will be considered in the fu- 
ture. 
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